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Abstract 

At a continuous transition into a nonunique absorbing state, particle systems may 
exhibit nonuniversal critical behavior, in apparent violation of hyperscaling. We pro- 
pose a generahzed scahng theory for dynamic critical behavior at a transition into an 
absorbing state, which is capable of describing exponents which vary according to the 
initial configuration. The resulting hyperscaling relation is supported by simulations 
of two lattice models. 
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Nonequilibrium phase transitions con- 
tinue to elicit great interest in the phys- 
ical and biological sciences. In hopes of 
better understanding nonequilibrium crit- 
ical phenomena, models exhibiting phase 
transitions into an absorbing state are under 
intensive study in statistical physics. The 
dichotomy between an absorbing (dead, in- 
active) state and an active one arises nat- 
urally in such diverse areas as catalysis 
epidemiology 0, ^, and the transition to 
turbulence 

The essential features of the phase tran- 
sition are typified by the contact process 
(CP) 0. In the CP, each site of the lat- 
tice Z'^ is either occupied or vacant. Occu- 
pied sites become vacant at unit rate, whilst 
a vacant site i becomes occupied at rate 
Xqi, with Qi the fraction of occupied near- 
est neighbors of i. Evidently, the vacuum 
is absorbing. The growth rate A determines 
the ultimate survival of the population: for 
A < Ac the vacuum is the unique steady 
state, but for A > Ac(— 3.298 in one dimen- 
sion), there is also an active stationary state 
characterized by a nonzero particle density 
p oc (A — Ac)^. The transition at Ac is a 
nonequilibrium critical point, belonging to 
the universality class of directed percolation 
(DP) i, |, 0] and Reggeon field theory [§. 
Indeed, studies of a host of models provide 
ample support for the conjecture that con- 
tinuous transitions into a unique absorbing 
state generically fall in the DP class 0, |10|. 



fall in the DP class 
behavior is concerned 



as far as static critical 
g, ig. (In all of 



The situation regarding models with mul- 
tiple absorbing states is more complex. On 
one hand, studies of some two-dimensional 
catalysis models yield critical exponents dif- 
ferent from those of DP [jll], . On the 
other, the one-dimensional pair contact pro- 
cess (PCP) and dimer reaction (DR) clearly 



these models, the number of absorbing con- 
figurations grows exponentially with lattice 
size.) The dynamic critical properties of the 
PCP and DR were found, surprisingly, to 
be nonuniversal, the exponents depending 
upon the nature of the initial configuration 
T5[| . While this variation is quite regular, it 



appears to violate a basic hyperscaling re- 
lation amongst the exponents 6, rj, and z 
(defined below), suggesting a breakdown of 
the well-established scaling theory. 

This apparent breakdown has prompted 
us to reexamine the scaling hypothesis for 
models with multiple absorbing states. We 
arrive at a scaling theory in which addi- 
tional exponents are expected to depend 
upon the starting configuration, and in 
which the exponents satisfy a generalized 
hyperscaling relation. The latter is verified 
in simulations of the DR and of a new model 
called the threshold transfer process. In 
what follows we define the models, present 
the scaling theory, and report the numerical 
evidence supporting it. 

Simulations of the threshold transfer pro- 
cess (TTP) permit us to study nonuniversal 
critical spreading in the context of a three- 
state model, providing a check on the ro- 
bustness of earlier findings [15]. In the TTP 
each site may be vacant, or singly or dou- 
bly occupied, corresponding to o"j = 0, 1 or 
2. In each cycle of the evolution, a site 
i is chosen at random. If at = 0, then 
(Tj — s> 1 with probability r; if (Xj = 1, then 



(Tj — > with probability 1 



r. 



'0" and 



"1" sites are left unchanged with probabil- 
ities 1 — r and r, respectively.) In the ab- 
sence of doubly-occupied sites, we have a 
trivial dynamics in which a fraction r of the 



1 



sites have cTj = 1, in the steady state. How- 
ever, if cTj = 2, the particles may move to 
neighboring sites. If (Tj_i < 2, one particle 
moves to that site; independently, a particle 
moves from i to i + 1 if a^+i < 2. cij is di- 
minished accordingly in this deterministic, 
particle-conserving transfer. Survival of the 
doubly-occupied sites hinges on the process 
(1, 2, 1) — i> (2, 0, 2) (their number decreases 
or remains the same in all other events), 
and so depends upon the parameter r, which 
controls the particle density. The set of con- 
figurations devoid of doubly-occupied sites 
comprises an absorbing subspace of the dy- 
namics, which can be avoided only if r is suf- 
ficiently large. Thus we identify the density 
of doubly-occupied sites, p2, as the order pa- 
rameter of the threshold transfer process. 

We note in passing that the TTP bears 
some resemblance to a sandpile model de- 
vised by Manna [0, and to a forest fire 
model (FFM) proposed by Bak, Tang and 
Weisenfeld fl^, Under the correspon- 
dence: 2 burning tree, 1 ^ live tree, and 
<r^ ashes, the process 1,2, 1 — 2,0,2 de- 
scribes a burning tree setting its neighbors 
on fire, and r represents the rate at which 
new trees emerge from the ashes. But the 
TTP permits doubly occupied sites to arise 
only via transfer; there is no "lightning" 
process, as in the FFM. A further contrast 
is that a rule such as 0, 2, — 1, 0, 1 has no 
analog in the FFM. Despite certain common 
aspects, our model is therefore very different 
from the FFM. 

The dimer reaction (DR), introduced in 
Ref. 
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is a lattice model in which sites 
may be either vacant or (singly) occupied; 
particles may not occupy adjacent sites. In 
each step of the process a site i is selected at 
random. If i is occupied, or is blocked by a 



neighboring particle then nothing happens. 
But if i is open (i.e., i and its neighbors are 
vacant), a new particle appears, which may 
remain at or react with another particle, 
depending on the occupancy of the nearby 
sites: 

(i) If at least one of the second neighbors, 
i — 2, and i + 2, are occupied, then with 
probability 1—p there is a reaction between 
the new particle and its neighbor (chosen at 
random if there is a choice), which removes 
them both; with probability p there is no 
reaction and the new particle remains. 

(ii) If both second neighbors are vacant, 
but at least one of the third neighbors 
is occupied, then a reaction with a third- 
neighbor particle may occur, with probabil- 
ities as in case (i). 

{Hi) If none of the second or third neigh- 
bors are occupied, the new particle remains 
at site i. 

In the DR, any configuration devoid of open 
sites is absorbing. The order parameter is 
the stationary open site fraction, p^. 

Both the TTP and the DR exhibit con- 
tinuous phase transitions to an absorbing 
state marked by a vanishing order param- 
eter at a critical value of r or of p. In this 
work we are concerned with critical spread- 
ing, that is, the evolution of a critical system 
from a nearly-absorbing initial configura- 
tion. The exponents describing this spread- 
ing are nonuniversal, i.e., they depend upon 
the particle density in the initial state . 
Before reporting our numerical results, we 
present a scaling theory for such processes. 

Following Grassberger and de la Torre 
I Toll , we consider an ensemble of trials, all 
starting from the same initial configuration: 
a single seed in an otherwise absorbing con- 
figuration. (For the contact process, this 
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means one particle in an otherwise empty 
lattice; for the DR, one open site; for the 
TTP, one doubly-occupied site.) Let p{x,t) 
denote the local order-parameter density, A 
the distance from the critical point (A = 
A — Ac in the contact process). In the critical 
region the system is characterized by a cor- 
relation length oc A"'^-^, and relaxation 
time r oc A"*^!!. At the critical point the 
asymptotic evolution is described by power- 
laws; for A 7^ the power laws are modified 
by scaling functions which depend upon the 
dimensionless ratios x/C,± and t/r. Thus the 
survival probability - i.e., that a trial has 
evaded the absorbing state, p{x) = - is 
expected to follow 



so that P oc at the critical point. The 
order parameter density (averaged over all 
trials) is 



(2) 



where the x-dependence reflects symmetry 
and power-law critical spreading from the 
seed at x = 0. For A = 0, one finds upon 
integrating Eq (0) over space, that the mean 
population n(t) oc f, whilst the second mo- 
ment implies that the mean-square spread 
of the population i?^(t) = {x'^)t oct^. The 
exponents S, t] and z characterize critical 
spreading; several relations connect them 
with other exponents. 

Consider first the CP, for which the 
ultimate survival probability. Poo = 
limt^oo P(t) = p, the stationary particle 
density |TP[. Existence of the limit requires 
0(x) OC x^'^^ii for large x, and p oc A^ then 
implies the scaling relation 



For A < and (large) fixed t we ex- 
pect p{x,t) ~ e~^/^, which implies (since 
^ oc A-"^), that for v < 0, F{u,v) oc 
exp{— const. y/u\v\'^^) . In order that ^ be 
time-independent, we must have 



(4) 



Finally, note that for A > 0, the local den- 
sity at any fixed x, in a surviving trial, must 
approach p as t — > cxd. Since p{x,t) repre- 
sents an average over all trials, we have 



p(x,t)^PooA^ocA2/^, 



(5) 



as t ^ oo, which implies that P(0, v) oc f 
for large v. Existence of a stationary state 
then requires that the exponents satisfy the 
hyperscaling relation 



A6 + 2r] = dz. 



(6) 



We turn now to models such as the pair 
contact process (PCP) |T^, the dimer reac- 



6 = (3/u\\. 



(3) 



tion and the threshold transfer process, 
which possess a multitude of absorbing con- 
figurations. Absorbing configurations in the 
PCP and DR can have various particle den- 
sities; the analogous variable in the TTP is 
the density pi of singly occupied sites. We 
refer to this aspect of the (near-absorbing) 
initial configuration in a critical spreading 
process as the "initial density," 0j. One 
value of (pi is special in these models: the 
"natural" particle density 0c of the qua- 
sistationary critical process. (For the DR, 
0c ^ 0.418; for the TTP, 0c ^ 0.69.) Sim- 
ulations indicate that in each of these one- 
dimensional models, the static critical be- 
havior belongs to the directed percolation 
class, but that the exponents 6, rj and z are 
nonuniversal, varying continuously with ini- 
tial density. (The critical point Ac, by con- 
trast, does not change as 0i is varied.) Only 
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when 0j = 0c do the exponents assume DP 
values, and only then do they satisfy Eq d). 
Rather than interpreting this violation 
of hyperscaling, we shall argue that in these 
models the scaling hypothesis must be mod- 
ified, leading to a generalized hyperscaling 
relation. 

We assume that as in models with a 
unique absorbing state, the order parame- 
ter density has the scaling form 

p{x,t) ~ t^'-''''/^G{xyt'',At'^'''\), (7) 

where the primed exponents are functions of 
(pi. Similarly, we suppose the survival prob- 
ability follows 



(8) 



Since the stationary distribution is unique, 
we have as before that 



(9) 



as t ^ oo, with /5 the usual directed perco- 
lation exponent. However there is no reason 
to suppose that Poo(0i) oc p, when 0j ^ 0c- 
In fact if this were so, we would have 5V|| = 
f3, implying that the primed exponents sat- 
isfy Eq(|D. Since they do not, we conclude 
that the exponent governing the ultimate 
survival probability must also depend upon 
0i, i.e.. Poo oc A^', with /?' = (5V[|. By the 
same arguments applied to the CP we find 



2^1/^11, 



(10) 



where we have introduced exponents z/'j and 
z/^ which govern the mean lifetime and spa- 
tial extent of a cluster grown from a single 
seed. The asymptotic behavior of the order 
parameter density is now p{x, t) — >■ A^~^^ , 



and G{0,y) oc ?/^+^' for large y, which im- 
plies the generalized hyperscaling relation 



2{l + ^)S' + 2v' 



dz'. 



(11) 



We have verified Eq (0) in simulations of 
the dimer reaction and the threshold trans- 
fer process. Using time-dependent simula- 
tions (for t < 2000, and samples of 10^ to 
8 X 10^ trials), we determined the critical 
point of the TTP as = 0.6894(3). Anal- 
ysis of steady-state data for pg, as shown in 
figure 1, then yields P = 0.279(5), in good 
agreement with the value for DP in 1 -f 1 
dimensions, (3 = 0.2769(2) |0|, ||]. The 



exponents 5', rj' and z' may be determined 
from simulations at Vc, using an initial con- 
figuration very close to the absorbing state. 
We studied various initial densities, includ- 
ing 0j = 0.69, the natural value. The ul- 
timate survival probability exponent jS' was 
determined from similar studies, using r val- 
ues slightly above critical. The simulations 
begin with one doubly-occupied site at the 
origin; the remaining sites are taken as oc- 
cupied or vacant, independently, with prob- 
abilities 0j and 1 — 0j, respectively. The 
dynamics is restricted to an active region 
defined as follows. Let Aj be the set of all 
sites which are doubly-occupied or have a 
doubly-occupied neighbor, after the i^^ step 
of the trial. (Aq comprises the origin and its 
neighbors.) The site to be updated at step 
2 + 1 is selected at random from U*^o ^i- 
Thus the evolution proceeds on an expand- 
ing set within the "light-cone" emanating 
from the origin. As in the DR and the PCP, 
distant sites are not updated until the active 
region reaches their neighborhood. Figure 2 
shows a local-slope analysis for 6, i.e., a plot 
of 6{t) = \n[P{mt)/P{t)]/\nm vs. r\ for 
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various initial densities. (In this study we 
used m = 5.) In figure 3 we show typical 
results for the ultimate survival probability, 
leading to an estimate for 

The simulation procedure for the DR is 
described in Ref. [^. On the basis of 
more extensive studies of the half-life r, on 
lattices of up to 1000 sites, we now find 
Pc = 0.26401(2), consistent with the ear- 
lier result of 0.26400(5). We used p values 
slightly below critical {pc — P ^ 0.05) in de- 
terminations of P' at the four initial densi- 

According to Eq R 



ties studied in Ref. [15 



in a plot of P = A'^'P vs. t = M '^w, data 
for various A (for a particular c/jj), should 
fall on a single curve. Figure 4 shows a rea- 
sonably good collapse of the data for each 
of four initial densities. 

Our results for the exponents in the TTP 
and DR are given in Table I, together with 
a test of the new hyperscaling relation, Eq 
(pT]). Evidently it is confirmed to within 
the precision of the data. (By contrast, the 
DP hyperscaling relation, Eq (^, is clearly 
violated.) Thus the spreading and survival 
exponents for transitions into a nonunique 
absorbing state may be described using the 
conventional sort of scaling theory, properly 
generalized to allow for a dependence upon 
the initial density. 

The dependence of 5', rj' and /?' upon ini- 
tial density is quite pronounced, that of z 
much weaker. We have made no determi- 
nation of z/^, and our results for z/'j, which 
come solely from the relation z/|| = 
show no significant variation with (pi. (We 
find v\\ = 1.80(6) and 1.76(6) for the TTP 
and the DR, respectively; the DP value is 
1.74(1).) In light of Eq (|TOD, it appears that 
i/_L is not strongly dependent upon initial 
density either. 



Further examination of the data indicates 
that S'+t]' is also very nearly constant. This 
is clear from the plot rj' vs. 5' for all three 
models (TTP, DR, and PCP) shown in fig- 
ure 5. (The slope of the linear best-fit is - 
0.995.) It is also worth noting that the expo- 
nents of the (two-dimensional) dimer-trimer 
model |]12| differ from those of DP, but that 
5 + Tj is again the same as in DP. (Sim- 
ulations of the dimer-trimer model yield 
5 = 0.40(1), T] = 0.28(1), compared with 
0.460(6) and 0.214(8), respectively, for two- 
dimensional DP [^.) Now S' + r)' is the ex- 
ponent governing the population growth in 
surviving critical trials. Its independence of 
(pi suggests that the asymptotic properties 
of a surviving trial are not affected by the 
initial density. This conclusion is strongly 
supported by the absence of any detectable 
shift in the critical point as (pi is varied. 
As further support, one may note that as 
t — i> oo, only a negligible fraction of a sur- 
viving cluster is actually in contact with the 
external density 0j. Deep inside the cluster, 
the particle density must approach the nat- 
ural value (pc- This point of view also implies 
that z', which describes surviving trials ex- 
clusively, should be constant. In fact, if z' 
and 6' + rj' are constant, then Eqs (|10|) and 
( |Tll) require that z/|| and z/^ are as well. We 
are led, by this line of argument, to a more 
economical description of critical spreading, 
in which all of the exponent variation fol- 
lows from a single cause: the dependence 
of the survival probability upon initial den- 
sity. The ensuing predictions regarding ex- 
ponents are consistent with our numerical 
results, except for a small variation in z' 
with (pi. One may argue, however, that the 
results for z' are affected by 0j-dependent 
corrections to scaling, and that a more pre- 
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cise numerical test is needed. In summary, 
we believe that the most natural and par- 
simonious interpretation is that the initial 
density influences the survival probability, 
but not the scaling properties of surviving 
events. 



The pair contact process, dimer reaction, 
and threshold transfer process all involve a 
second variable, 0, dynamically coupled to 
the order parameter. A quantitative theory 
of the dependence survival probability, and 
the associated exponents 5' and /?', upon the 
initial density (pi has yet to be devised. But 
we can offer some intuitive basis for under- 
standing nonuniversality by suggesting that 
in these processes, the initial density plays 
a role analogous to that of a marginal pa- 
rameter. Such parameters, invariant under 
renormalization group transformations, of- 
ten give rise to exponents which vary con- 
tinuously along a line of fixed points. In the 
present case, (pi represents a property of the 
medium into which the process grows, and 
which is never forgotten, since to survive, a 
critical process must repeatedly invade new 
territory. A RG transformation generally 
involves coarse graining (which conserves 
density in the large), and rescaling. Such 
a transformation may be expected to leave 
(pi, the density outside the active region, in- 
variant, while driving the correlation length 
of this region to zero. Indeed, the spreading 
exponents are insensitive to (short range) 
correlations in the initial state |jl5| . A more 
detailed understanding of nonuniversality in 
these models may emerge when a suitable 
renormalization group scheme is devised. 
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Figure Captions 



Figure 1. Steady state order parameter den- 
sity, P2, versus distance from critical point 
in the TTP. (Main graph: log-log plot; in- 
set: linear plot.) 

Figure 2. Local slope analysis of the sur- 
vival probability data for various 0j values 
in the TTP. 5 is estimated from the i — > 00 
intercept. 

Figure 3. Ultimate survival probability ver- 
sus distance from the critical point in the 
TTP, for initial density 0^ = 0.4. 

Figure 4. Scaling plot of the survival prob- 
ability in the DR. -F: A = 0.05; x : 

A = 0.02; □: A = 0.01; o: A = 0.005; o: 
A = 0.002. Initial densities (top to bottom) 
(t)i = 1/2,0.418,0.38, and 1/3. 

Figure 5. rj' vs S' for the PCP, DR, and 
TTP. The slope of the best-fit straight line 
is -0.995. 
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Table 1: Critical exponents of the thresh- 
old transfer process and the dimer reaction. 
Figures in parentheses denote uncertainties 
in the last figure (s). 



0i 


S' 


1]' 


z'/2 




rj'+{l + ^)6'-z'/2 


Threshold Transfer Process 


0.75 


0.136(1) 


0.347(4) 


0.632(7) 


0.250(5) 


0.00(1) 


0.69 


0.161(2) 


0.319(3) 


0.632(7) 


0.279(5) 


0.00(1) 


0.60 


0.192(2) 


0.288(3) 


0.630(7) 


0.356(5) 


0.00(1) 


0.50 


0.227(2) 


0.246(2) 


0.623(7) 


0.426(5) 


0.00(1) 


0.40 


0.270(3) 


0.204(2) 


0.622(7) 


0.497(5) 


0.00(1) 


0.31 


0.299(3) 


0.169(2) 


0.617(7) 


0.556(6) 


0.00(1) 


0.30 


0.303(3) 


0.170(2) 


0.621(7) 


0.567(6) 


0.00(1) 


0.25 


0.316(3) 


0.161(2) 


0.624(7) 


0.591(6) 


0.00(1) 


0.20 


0.342(3) 


0.133(1) 


0.622(7) 


0.640(6) 


0.00(1) 


0.10 


0.371(4) 


0.097(1) 


0.615(7) 


0.705(7) 


0.00(1) 


Dimer Reaction 


0.333 


0.107(2) 


0.362(3) 


0.634(3) 


0.182(10) 


0.00(1) 


0.380 


0.133(2) 


0.327(3) 


0.629(5) 


0.241(6) 


-0.02(1) 


0.418 


0.158(2) 


0.302(4) 


0.626(3) 


0.275(2) 


-0.01(1) 


0.500 


0.205(5) 


0.250(5) 


0.620(3) 


0.357(10) 


-0.01(1) 
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